Resonances, instabilities, and structure selection of driven Josephson lattice 

in layered superconductors 
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We investigate dynamics of Josephson vortex lattice in layered high Tc superconductors at high 
magnetic fields. It is shown that the average electric current depends on the lattice structure and 
is resonantly enhanced when the Josephson frequency matches the frequency of the plasma mode. 
We find the stability regions of moving lattice. It is shown that a specific lattice structure at given 
velocity is uniquely selected by the boundary conditions: at small velocities periodic triangular 
lattice is stable and looses its stability at some critical velocity. At even higher velocities a structure 
close to a rectangular lattice is restored. 



Transport properties of layered superconductors, such 
as Bi2Sr2CaCu202; (BSCCO), in magnetic field parallel 
the layers are determined by the dynamics of Joseph- 
son vortex lattice (JVL). ||l|-g[ Moving JVL generates 
travelling electromagnetic wave in the medium. Similar 
to the Eck resonance in a single junction 0], one ex- 
pects a strong resonance emission when the velocity of 
the lattice matches the plasma wave velocity. ^^ In the 
current-voltage (I-V) dependence this resonance is seen 
as a strong enhancement of current at fixed voltage. In 
contrast to a single junction, JVL in layered supercon- 
ductors has soft degrees of freedom related to the phase 
shifts between different layers. This leads to a rich variety 
of dynamic states observed in numerical simulations. R] 
In particular, a periodic lattice corresponds to a constant 
phase shift between neighboring layers, ranging from for 
the rectangular lattice to n for the static triangular lat- 
tice. Moving lattice generates electromagnetic wave with 
the c-axis wave vector selected by the lattice structure. 
Since the velocity of a plasma wave propagating along 
the layers depends on this wave vector |g|, the resonant 
velocity of the lattice depends on its structure. 

In this Letter we investigate the stability of moving 
JVLs and the evolution of structure as a function of its 
velocity for large size samples. [|| We show that a specific 
lattice structure at given velocity is uniquely selected by 
the boundary conditions. At small velocities a periodic 
lattice is stable. The phase shift between neighboring 
layers smoothly decreases with the increase of velocity, 
starting from tt for a static lattice. At some critical ve- 
locity the lattice becomes unstable. At even higher veloc- 
ities a periodic lattice with the phase shift smaller than 
tt/2 is restored again. 

Consider a layered superconductor in a magnetic field 
applied along the layers (y axis) with transport current 
flowing across the layers (z axis) . The dynamics of such 
a system is described by the inductively coupled sine- 
Gordon equations for the gauge-invariant phase differ- 
ences On- pO| | Introducing the reduced time r — cupt and 
in-plane coordinate u = x/{js), we express currents jx.n, 
jz,n and electric fields E^.n via 0„ and the reduced mag- 
netic fields hn = i?n27r7A^j/$o as 



jz,n = jj (sm6'„ + Vabd9n/dT) , E^^n ~ EpdOn/dr, 

_. _ ■ (,^ d_\ (de^_K 

Jx,n+1 Jx,n - Jab \i^ + l^ab g^ j \ Q^ ^ 

where ujp is the Josephson plasma frequency, 7 is the 
anisotropy factor, s is the interlayer spacing, jj is the 
Josephson current, Ep = ^quJp/(2ttcs), I — Xab/s, and 
jab = c$o/(87rA^f,7s). Vc = -iirac/iecUip) and Vab = 
Airaab^ab^p/'^'^ ^^^ the dissipation parameters, which are 
determined by the quasiparticle conductivities ac and aab 
and are connected as Uc/vab = Oc'i^ joab- Using the above 
relations we derive from the Maxwell equations the cou- 
pled equations for Q^ and /i„ 
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where V^/i„ = /i„+i -|-/i„_i — 2/i„. Equivalent forms of 
these equations have been derived in Refs. ||ll|,p|. In the 
case of negligible in-plane dissipation (vab = 0) they can 
be reduced to equations containing only 9„. ||l^ Taking 
parameters typical for underdoped BSCCO at T«50K: 
7 = 500, Aah = 240nm, s = 15A, cra6 = 2-104 [rj-cm]-i, and 
CTc = 2 • 10~^ [n ■ cm]~^, one obtains estimates: z/a6~0.1 
and i^c ~ 0.002, which we will use in numerical compu- 
tations. For these parameters the dissipation in a wide 
range of electric and magnetic fields is mainly determined 
by the in-plane channel. 

We consider the situation when all junctions are filled 
by JVs. This is always the case at high magnetic field B 
above ^o/^-^ls"^ [Q '^^ ^^ high currents, when all junc- 
tion are driven into the resistive state. Neglecting self- 
field of the transport current, we represent a solution to 
Eqs. (|l|) and (||) for the resistive state in the form 

Onir, u) = loet + knu + (/>„(r) + Onir, u) (3) 

hn{T,u) = h + hn{T,u). (4) 

Here the Josephson frequency ue is determined by the 
electric field, uje = E^/Ep, and the wave vector ku is 



determined by the magnetic field, kH = 2TrHjs'^/^o = 
h/P. 0n{T,u) and hn{T,u) are the oscillating phase and 
magnetic field induced by Josephson coupling, 0„,hn oc 
sm{uJET + knu + a„). Solution (|^) corresponds to the 
lattice moving with velocity (jjE/ku- The structure of 
the lattice is determined by the phase shifts (/)„. Without 
Josephson coupling the system is degenerate with respect 
to arbitrary phase shifts 0„. Josephson coupling elimi- 
nates this degeneracy. It either leads to slow dynamics 
of phase shifts 0„ or selects a certain steady state struc- 
ture. The equations for </)„(r) can be obtained using an 
expansion with respect to Josephson coupling and aver- 
aging over fast degrees of freedom. This procedure will 
be described in detail elsewhere. |l3|] 

For periodic JVL one has 0„ = nn. In such a state 
the chains of JVs in the neighboring layers are shifted 
by the fraction k/2'k of the lattice constant (see inset in 
Fig. H). In particular, k = corresponds to a rectangu- 
lar lattice, and k — t: corresponds to a triangular lattice 
(ground state at ue — Q)- In first order with respect to 
Josephson current, the oscillating phase 6'„ is given by 
9n = Im [Q{k) exp (i{ijjET + k^u + ktt,))] with 



g{K). 



UJe — IVcUOe 



kjj (1 +iVahtOE) 



2(1-COSk) + {l+iVabt^E)/l'^ 



To second order we obtain the average reduced Joseph- 
son current ij = ij{K, kH,(^E) — (sin0„(r, u)). |0| In 
the case of not very large in-plane dissipation Vab^E ^ 
2P(1 — cosk) + 1, it has the resonant dependence on the 
Josephson frequency loe (electric field) 
I 



iOEv{n) 



where Wp(k) = ku (2 (1 — cos k) + p) is the plasma 

frequency at the wave vector k and 



v{k) = i^c + 



2(1 - cos k) kfjl^ab 

(2(1-cosk) + 1//2)2 



is the dissipation parameter of the plasma mode. lO] 
When the frequency loe matches the corresponding 
plasma wave frequency (^^(k), a resonance enhancement 
of the current is expected. The moving lattice generates a 
travelling electromagnetic wave with the Poynting vector 
P. For the in-plane component of P we obtain 



Prr. 



PabUJEUJpiK)/kH 



{uj%-ul{n)Y + {iUEu{K)y 



:, Pab 



32tt^X^sj' 



(5) 



For typical BSCCO parameters at low temperatures 
(A = 200 nm, 7 = 500, Wp/27r = 150 GHz) we obtain 
Pab «125W/cm2. 

To investigate the stability of JVL we perturb the 
lattice solution as (pniT) — nn + I{.e[uq exp{a{q)T + 
iqn)] and derive from Eqs. (|l],||) the eigenvalue a{q) = 
a{q,K,kH,ujE) 

I , 1 g{K + q) + g*{K-q)-2^e[g{K)] 



4 ^c + i [J^ (k + g, k) - T* (-K -f- q, K)] 



with 



^{q,qi) = -2LJi 



ll^Jjl^ab 



2(1 — cosqi) 



2(1 — cos (7) 



W;;""^ 2(1^ cos gi) 



l + iVgb^E I 
P / 



g{q)g{qi)- 



The lattice is stable if there is no exponentially growing 
solution in the whole q-interval, i.e., Re[a((7)] < 0, for 
< <7 < TT. The onset of instability is characterized by 
the wave vector qi of the most unstable mode. There are 
three special cases: the long-wave instability qi = 0, the 
short-wave instability qi = tt, and the instability with 
< qi < TT. One can also distinguish absolute and con- 
vective instabilities. Pq] The spectrum a{q) has an im- 
portant symmetry property a(q = n, k = n/2) = 0. This 
means that a stable region cannot cross the line k ~ -k 12 
and it is impossible to evolve continuously from the static 
triangular lattice to the fastly moving rectangular lattice 
without intersecting an instability boundary. The trian- 
gular lattice always looses its stability before reaching the 
resonance uje — fcff/2, at the frequency WA(fc_ff) which is 
determined by a simple analytical equation 



= -UJA V, 



^H'^ab 



1 + J^ah^A ^ l^abi^A 



To explore the stability of JVL we calculated numeri- 
cally Re[a(q)] (ph throughout the loe — h plane and found 



the lines at which either Re[a(g)] = at finite q or 
Ke[cPa{q)/dq'^] = at q = 0. The stability diagram 
for representative parameters i>c — 0.002, Vab — 0.1, and 
kn = 8 is shown in Fig. P. For these parameters we 
found three stability regions at moderate loe'- (i)the low- 
velocity region located below the resonance line and at 
tt/2 < k < tt, (ii)thc high-velocity region located along 
the resonance line with k approaching with increase of 
lue^ and (iii)the region located above the resonance line 
and at 7r/2 < k < tt (this region disappears at higher 
fields). At the boundary of the first region the lattice ex- 
periences a long-wave instability for k > 2.04. At smaller 
K the instability occurs at finite wave vector q — qi and 
qi grows continuously with decrease of k. 

Selection by the system of a specific wave number k 
from the continuous spectrum is not directly related to 
stability of the corresponding periodic structure. For the 
static case the structure is selected by the minimum en- 
ergy condition. Such a condition is absent in the dynamic 



case. In this case a specific JVL structure can be deter- 
mined by the boundary. To demonstrate this, we con- 
sider a semi-infinite stack of junctions with n = 1,2... 
separated by a sharp boundary from the medium with ar- 
bitrary electromagnetic properties. To derive equations 
for the phase shifts </)„ for such system we have to find a 
solution of the linear equations without Josephson cou- 
pling taking into account the boundary conditions. For 
plasma wave with given frequency uj — uje and wave vec- 
tor along the layers k = kn the oscillating phases (0„) 
and magnetic fields (/i„) in the junctions can be written 
as 

On, hn oc exp(— i(7+ri) + B e-xp(iq^n) , at n > 1 (7) 

where g^ = q-f-{k,Lu) is (complex) wave vector given by 



cos (7+ = 1 



/C^ (1 -I- Wabi^) . 1 + il^ab^ 



2(^2 



IVcl-^) 



2/2 



(8) 



with Im[g+] > 0. The properties of the boundary are 
completely characterized by the complex amplitude of 
reflected wave B = B{k,uj), which has to be found by 
matching the solution ^ with electromagnetic oscilla- 
tions in the medium at z < 0. In general, B{k,Lu) can be 
a complex number with arbitrary absolute value. Only in 
the simplest case of vanishing dissipation and propagat- 
ing wave (lm((7+) = 0), B{k, to) determines a conventional 
reflection coefficient, R{k,uj) = \B{k,uj)\ , and has prop- 
erty \B{k,uj)\ < 1. A large class of boundaries, including 
boundary with free space, is well described by the ideal 
reflection B = —1. Averaging with respect to oscillating 
phases and field we derive the following equation for the 
steady state phase shifts </>„ 



1 °° 

- ^ Im [G{n, m) exp (- 



</'m))] = «J, 



where 



G{n, m) = Go{n — m) + BGo{n + m), 
Goin) = GQ{n;kH,ujE) = 



(9) 



(10) 
(11) 



Sn 



Je - il^ci^E 



kjj (1 + JUgbLOE) expig+l? 
{uolj^-iVcOJEf 2ising4 



The second term in Eq. (|0|) describes the surface con- 
tribution and vanishes at n^m —f oo. The solution, 
corresponding to the lattice in the bulk, has the form 
(j)n — Kn + Un, where u„ is the surface deformation, 
M„ ^ at n — > oo. Substituting this ansatz into Eq. 
(0) we obtain a nonlinear degenerate system of equations 
for Un- The solution for u„ exists only for special values 
of K, i.e., the bulk structure is selected by the boundary 
conditions. For such selected states z component of the 
Poynting vector Pz is always directed from the boundary 
towards the bulk of the sample. pM A similar pattern se- 
lection mechanism is relevant for various nonequilibrium 



systems. [|^ At small uje the system becomes linear with 
respect to m„ and its solvability condition can be found 
analytically. In the case B ^ — 1 we obtain deformation of 
the lattice at small velocities k ~ n — (vab + Ivc/kf^) lue- 
In finite system with identical boundaries the configura- 
tion is typically symmetric because each boundary selects 
the same wave number. The waves collide in the middle 
forming a phase defect, jl^ 

To find the steady state configurations at all velocities 
we solved Eq. (0) numerically taking B = —l and using the 
same parameters (i^c = 0.002, i/at = 0.1). The dependence 
k(uje) obtained from these solutions for kH = S is shown 
in Fig.|l| together with the stability regions. At small ve- 
locities the lattice experiences smooth evolution of struc- 
ture with the lattice wave vector k decreasing from n at 
zero velocity to smaller values until it hits the instabil- 
ity boundary. At higher velocities the stable lattice with 
K < 7t/2 is restored. The structure continues to evolve 
smoothly towards the rectangular configurations with in- 
crease of velocity. Near the line k — tt/2 we observe tran- 
sition to the double-periodic lattice 0„ —TTn/2 + {-'l)"v, 
which becomes stable at high fields. 

Fig. |2| shows evolution of the current-voltage depen- 
dence and the in-plane Poynting vector (^ with in- 
crease of magnetic field. Up to magnetic field H « 
\/oab/(j^<7c)^()/{njs^) {kn ~ 16), the current-voltage 
dependencies have two stable branches, corresponding 
to moving periodic lattices, separated by broad insta- 
bility region where periodic JVL can not exist. In this 
regime two lattice solutions exist within a finite range 
of currents. Such coexistence is facilitated by the high 
in-plane dissipation, which causes strong dependence of 
the JVL velocity on its structure. Within the coexistence 
region one can expect a family of intermediate states, in 
which the system is split into two (or more) domains mov- 
ing with different velocities separated by a phase defect 
(dynamic phase separation, see Fig. (^). The phase- 
separated states give the most natural interpretation of 
the multiple I-V branches observed by Hechtfischer et 
ai, who studied transport properties of JVL in BSCCO 
mesas at high magnetic fields, yj This interpretation can 
be verified by measuring the spectrum of microwave ir- 
radiation emitting by the stack. Instead of a single peak 
located at the Josephson frequency corresponding to the 
average voltage, the spectrum of irradiation should con- 
tain two peaks corresponding to the "fast" and "slow" 
states. 

In conlusion, we investigated stability and boundary 
structure selection of the driven JVL. We found two ma- 
jor stability regions, separated by unstable region: the 
low-velocity region correspondsto moving structure close 
to triangular lattice and the high-velocity region corre- 
sponds to an almost rectangular lattice. 
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FIG. 1. Stability regions of moving Josephson lattice 
in the plane wb-k calculated for representative parameters 
Vc ~ 0.002, Vab = 0.1, and kn ~ 8. Grey regions corre- 
spond to unstable lattices. Sections of the boundaries marked 
by black correspond to the long-wave instability. Sections 
marked by light grey correspond to the instability with q = n. 
Line starting at (0, n) shows dependence of the lattice wave 
vector on the frequency uje selected by the boundary with 
free space. We also show the resonance line, corresponding to 
matching between the Josephson frequency and the frequency 
of the plasma wave at the wave vector k. Inset sketches a 
steady state corresponding to a regular lattice. 
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FIG. 2. (a) Current-voltage characteristics at different 
magnetic fields, for kn = 6, 8, 12, and 16 (for 7 — 500 
H ~ kn ■ 0.3T). The dependencies are obtained using numer- 
ically computed steady states with parameters i/c = 0.002, 
Vab ~ 0.1. Thick lines show stable branches and thin lines 
show unstable branches. The branches, marked by dashes, 
correspond to the double-periodic lattices, (b) Electric field 
dependencies of the Poynting vector along the layers for elec- 
tromagnetic wave generated by the moving lattice (see Eq. 
(|)). The scale Pat is defined in Eq. (||). 
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FIG. 3. Multibranch structure of the current-voltage char- 
acteristic due to dynamic phase separation. Two states cor- 
responding to slow lattice motion (velocity Vs) and fast lat- 
tice motion (velocity Vf) coexist within the current range 
marked at the vertical axis. In this region the intermediate 
phase-separated states exist, in which the system is split into 
lastly and slowly moving regions. The intermediate branch 
corresponding to one of such states is shown by dotted line. 



